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Abstract - In this paper, natural circulation phenomena in a pressurized water reactor during a small break loss of coolant accident is 
investigated analytically. By using a one-dimensional analyses and quasi-steady state for two-phase flow, the values of the mass flow 
rate versus the primary coolant pressure for different core power are determined. It is found that the mass flow rate is reduced suddenly 
when the two-phase natural circulation is started. Also, it is found that the effect of the core power on mass flow rate is insignificant in 
the two-phase natural circulation. The analytical results of this paper are compared with previous experimental results and show a very 
good agreement. 
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1. Introduction 

It is strongly believed that the natural circulation in a PWR type system can provide an effective mechanism for the 
rejection of core decay heat to the secondary side of the steam generator. The physics of circulation process is a manometric 
balance between the hydrostatic head in the down-comer driving the core inlet flow, and the pressure losses incurred in 
venting the resulting two-phase flow from the core outlet to the heat exchanger [1, 2]. 

In the literature, there are different studies in natural circulation. Gruszczynski and Viskanta studied the flow and heat 
transfer for water in a rectangular loop in which tube bundles served as the source of heat in one vertical leg and as the heat 
sink in the other vertical leg. The purpose of the work was to gain improved understanding of the behaviour of a natural 
circulation loop containing tube bundles. Transients and steady-state tests were conducted and an unsteady one-dimensional 
formulation was used to model the flow. Detailed comparisons were made between the predicted and measured temperatures 
and good agreement was obtained under transient and steady-state conditions [3]. 

Steimke carried out an experimental and theoretical investigation of natural circulation of water in a loop where the flow 
channels are the annuli between three concentric cylindrical surfaces. An analysis was made based on a one-dimensional 
model. Results were obtained from the numerical solution to the conservation equation of mass, momentum and energy for 
steady laminar flow. Agreement between the calculated and measured temperatures, velocities, and Nusselt numbers was 
reasonably good [4]. Britt and Wood carried out experiments with several systems to determine the operational characteristics 
of a multiple-loop container used to cool nuclear reactor assemblies at the Savannah River plant [5]. Cheng and Rovang 
studied the steady-state operation and heat transfer characteristics of a closed loop, two-phase thermosyphon system for low 
grade heat recovery from liquid heat sources [6]. Burchill and Abramson [7], and, Duffey and Sursock [8] identified three 
distinct modes of natural circulation relevant to small breaks and transients for PWRs as a function of quality: 1) single- 
phase, 2) two-phase, and 3) reflux condensation. Abdulrahman has investigated the analytical solution of the energy equation 
in two dimensional steady state condition. He has determined the temperature distributions of the packed bed reactor in terms 
of the radial and axial dimensions, as well as the temperature distribution as a function of time. [9, 10]. 

The present work is an analytical investigation of a two-phase natural circulation system relevant to a PWR during a 
small break loss of coolant accident. The model of this paper utilizes the one-dimensional analyses and the quasi-steady state 
and solves analytically the conservation equations of mass, momentum and energy. In this paper, an expression for the two- 
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phase mass flow rate as a function of a primary pressure is developed, as well as the effect of core power on the mas 
flow rate is investigated. 


2. Steady State One-Dimensional Natural Circulation Flow 
For a flow in a one-dimensional positive s-direction, the steady state equations of continuity, momentum, and energy 
are respectively [11]; 


du udp 
== (1) 
ds pds 
du dP 7 (6) (2) 
er a ds s pi cos 
q 
habeas any 8 3 
pat ds ds 2 


where u is the flow velocity, p is the fluid density, P is the primary pressure, g is the gravitational acceleration, A is the 
cross sectional area of the flow channel, h is the specific enthalpy, and q is the heat added to (or rejected from) the coolant. 
In terms of the mass flow rate W, the flow velocity u is: 


u = — (4) 


For turbulent flow in a circular tube of diameter D, the viscous force per unit volume F; is; 


1 2 
p Se) (5) 
f= — A ^ 
D 
where the friction factor f is given by [12]; f = au Substituting Eqs. (1), (4) and (5) into Eq. (2) and integrate around 


the loop, the momentum Eq. will be; 


1 
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where v is the fluid specific volume and P is the primary pressure. The left hand side and the pressure term of Eq. (6) will 


vanish. For strongly turbulent flow, it can be assumed that the friction coefficient is constant for a given primary pressure. 
Noting that A and D are given by their average values over the loop and that dz = ds cos(@), then Eq. (6) becomes; 


ds — fog cos(@) ds (6) 


TAE 
_ _Z4A)W? fds (7) 
gf p dz = RD 5 


206-2 


Steam generator — 


A 


D 


HE 
H 


E HEHEHEH 


Core pX 
t 


La 


1 


Fig. 1: One-dimensional loop. 


From Fig. 1; 
Lz 0 -La 0 Lz Lz 
$ pdz= | pidz+| pidz+ | prdz+ | prdz =| paz- | p2 dz = (Pı — P2)Lz (8) 
0 Lz 0 -La 0 0 
and; 
ds 
po = bvds=ol, (9) 
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where p4 and pz are the fluid densities of the system up-flow and down-flow sides respectively, p4 and Pz are the average 
values of the fluid densities p4 and p, respectively, v is the average value of the fluid specific volume, L,, and Lg are the 
heights of the PWR and downcomer systems respectively, and L+ is the total circulation length. Substituting Eqs. (8) and (9) 
into Eq. (7) and rearranging, to get; 


5 (4f)W? 
A2 D 


2 2 = 10 
DL, = g(P2 — pi) Lz (10) 
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where D and A are the average diameter and cross sectional area of the flow channel respectively. Note that W has the 
dimensions of flow rate and fọ and Jọ are respectively the overall system average friction coefficient and density 
corresponding to the conditions at the beginning of NC. Substituting Eq. (11) into Eq. (10) and solving for the mass flow 
rate, to get; 
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(12) 


where f is the average friction factor. The value of pı is obtained from volume averaging of densities over the up-flow side 
of the PWR system elements, as follows: 


(13) 


where pc, Pu, Ps are the average densities of the core, sum of the hot leg and upper plenum, and steam generator respectively. 
Va, Vi» Vs, and V; are the volumes of the core, sum of the hot leg and upper plenum, and steam generator respectively and 


VY, =V4th,4+ s, The energy Eq. (10) can be written separately for each component of the loop. The term (- za) in Eq. 


(10) is equal to (-UmD,n(T — Tsec)) for the heat sink and zero for the (insulated) pipes. For the heated section, it depends 
upon the input power distribution. It is well known that the power distribution in a typical PWR core is far from being 
uniform. To first approximation, the distribution may be taken as sinusoidal. In this paper, to simplify derivations, uniform 
input power distribution is taken into consideration. For a uniformly distributed input power, the term (- za) is equal to 2o 
for the heated section. Hence, Eq. (10) can be written as: 


Qo/Le (heat source) (14a) 
Wz = < -Un Dsn (T; — Tsec) (heat sink) (14b) 
0 (insulated pipes) (14c) 


where Q, is the heat input power, L, is the length of the core, U is the average value of the overall heat transfer coefficient, 
D, is the steam generator diameter, n is the number of primary tubes of steam generator, Ts and Tec are the steam generator 
temperatures of the primary and secondary sides respectively. 


2.1. Two-Phase Loop Flow 

In natural circulation, the two-phase mode is regarded to begin when the bubbly flow is carried over into the down-flow 
side of steam generator and exit piping. Consider the loop to contain only two-phase flow, with unequal densities J4 and po, 
on the hot and cold sides, respectively (Fig. 2). 
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Fig. 2: Two-phase loop. 
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The solution of Eq. (14c) yields uniform enthalpy and hence quality in the hot leg pipe, which implies, that the outlet 
core and inlet steam generator qualities are at xg. Utilizing the definition of h = hr + xhgg in Eq. (14a), integrating, and 
using the boundary conditions (x = x; at z = 0 and x = x, at Z = Lo) the quality at the core outlet can be found as; 


Xo = Xi + (15) 


W hfg 
where x; is the quality at the core inlet, and hpg is the latent heat of vaporization. By solving Eq. (14b) at z = 5, the quality 


at the core inlet (steam generator outlet) can be found as; 


PE UnDsn(Tsat g Tec) Ls 


aT (16) 
W hrg 2 
Equations (15) and (16) are used now to determine the average density for each elemnet of a PWR loop. 
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where Pe, Py; Ps, P2 are the average densities of the core, sum of the hot leg and upper plenum, steam generator, down-flow 
side respectively, vp is the saturated liquid specific volume, and Ugg = Ug — vf. The average specific volume, v , in Eq. (12), 
can be obtained from volume averaging over the PWR system elements, thus; 


TV, + DyVy + Ds $ + TeV, 
V, + Vz 


P (21) 


where Vz, Pu, Vs, and V2 are the average specific volumes of the core, sum of the hot leg and upper plenum, steam generator, 
and down-flow side respectively. V., V}, and V;, are the volumes of the core, sum of the hot leg and upper plenum, and steam 


‘ V, V, ; SEEPS 2 ; 
generator respectively, V = V; + Va + re and Vz = = + Va. The magnitudes of ¥,, Du, Vs, and V2 can be determined from: 


Ufg Q 
A fg xo 
Uc = Vf tUg Xi aaa (22) 
2 W hrg 
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2W hfg 2 
U, = Uf + Ugg Xi (25) 
By solving equations (15) and (16), the following expression can be obtained: 
z- Ls 
Qo =U T Dsn (Tsat — Tsec) = (26) 
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Since, in nucleate boiling, U is a function of a primary pressure, P and the other parameters in equation (26) are constants, 
so that, it can be concluded that, for two-phase flow, the primary pressure, P, is constant and can be evaluated as the same 
value of the primary pressure at the peak mass flow rate. At this stage, one proceeds to evaluate the average densities in the 
hot and cold legs, p; and p, , and the average specific volume, v . The mass flow rate, W, can be obtained from equation 
(12). It should be kept in mind that the temperature difference across the core will be zero, since the inlet and outlet 
temperatures are at saturation temperature. 


2.2. Overall Heat Transfer Coefficient 
The overall heat transfer coefficient for the radial flow of heat across a circular tube based on the internal tube surface, 


Uj, is [13]; 


1 1 Ajtw Aj 
ee eS pe Ee (27) 
U; H i AmwKw AoHo 

where H; and H, are the heat transfer coefficient inside and outside the tube respectively, A; and A, are the internal and 
external tube surface area respectively, Amw is the mean wall area for heat transfer, and Ky, is the thermal conductivity of 
the wall tube. The heat transfer resistance by conduction can be neglected because of small thickness of steam generator 
tubes. For two phase flow in the steam generator, convective heat transfer on the inside of the tubes and nucleate boiling on 
the secondary side of the tubes will occur [14]. The heat transfer coefficients in the primary side of steam generator H;, have 
been taken, from forced flow correlations. From the above considerations, the average overall heat transfer coefficient for 


two-phase flow, U, can be written as; 
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= 1 A; 
U =—— + —*_ (28) 
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where; Hiz¢ is the heat transfer coefficient for two-phase in the primary side of steam generator for bubble flow, which can 
be taken from the following correlation [2]: 


Hing = 2.5551 x 10-3¢964461357 (T. at — Trec)? (29) 


Hozọ is the heat transfer coefficient for nucleate boiling of the secondary side of steam generator, which can be taken 
from the following correlation [2]: 


Hogg 2.5551 X 10° tet Ate tsee(T at — Teas) (30) 


r ; : : ; kW 
where Psec is the secondary side pressure in MPa, T in °C and H, in zog To evaluate the two phase Reynolds number, 


two phase dynamic viscosity, 2g, must be determined from the following equation: 


1 ee. ee 
— = +——, (31) 


oes ae : r . it 
where; X, is the average quality of stream generator which is equal to ately 


3. Results 


In this paper, in order to obtain the results of the NC parameters, iterations are used. The values of p, and f, for every 
reactor power Qo, can be determined when the primary pressure is at its initial value (P,=11.2 MPa) (see Table 1). 


Table 1: Values of , and f, as a function of core power. 


Q» (KW) 30 60 100 
Po 743.2258 734.6119 730.6987 


Ía 4.847 x 10° | 4.614 x 10° | 4.462 x 10° 


The results of this paper are compared with the experimental data from the Semiscale Mod-2A system for three different 
powers (30, 60 and 100 KW) [14]. The data required for the calculations in this paper that concern the Semiscale Mod-2A 
system are shown in Table 2 [2, 8, 14]. Figures 3 and 4 show the calculated mass flow rate of the natural circulation as a 


function of the primary pressure as well as the comparisons of the calculated results with the experimental results of 
Semiscale Mod-2A system. 


Table 2: Summary of the design and operating data for the Semiscale model. 


Parameter! Data Parameter Data Parameter; Data 
Ai 3.05 x 10 74 m? Ls 20 m Prec 6 MPa 
Ao 3.87 x 10-4 m? Ve 0.06 Tsec 275.64 °C 
Ds 0.0222 m Va 0.18 Tı Toutk 
tw 0.00124 m V; 0.22 Pı Poulk 
n 6 Va 0.29 W 1.35 Kg/sec 
Le 3.66 m full core power (scaled)2 MW 
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Fig. 3: Comparison of theory with Semiscale PWR data of Mass flow rate versus Primary pressure for different core powers (a) Q, = 
30 kW, (b) Q, = 60 kW, (c) Q, = 100 kW. 
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Fig. 4: Effect of core power on Mass flow rate versus Primary pressure. 


In PWR, NC starts when all the primary coolant is at single-phase. When the draining in the primary loop is continued, 
the lever of the coolant will decrease until it reaches the hot leg where combined single and two-phase flow in the NC will 
start. As more fluid is drained in the primary loop of the PWR, the bubbly flow will eventually be carried over into the steam 
generator down-flow side and exit piping. In this case, the whole primary coolant will be in a two-phase mode. 

In Fig. 3, it can be seen that in two-phase flow, the mass flow rate is sharply reduced as the overall driving head is 
drastically reduced by the presence of two-phase flow in the down-flow side of the steam generator. The cold leg temperature 
will remain at the same saturated temperature because the pressure will remain constant for two-phase flow. Since both the 
inlet and outlet temperatures are at saturation, then the core temperature difference will be reduced to zero. Figure 4 shows 
the effect of the reactor power on the NC parameters. In general, it can be seen that the overall trends for the three power 
levels are similar. The mass flow rate at the starting of the two-phase NC, is greater the higher the core power, while at the 
ending point is approximately the same. 

Table 3 shows the maximum error between the experimental and theoretical results of the primary pressure and the NC 
mass flow rate. The discrepancies in the results can be attributed to different reasons. The first one is the uncertainties in 
experimental data, which were estimated [8] as + 0.033 Kg/s for mass flow rate, and + 0.4 MPa for secondary pressure. The 
second reason is the model approximations. The one dimensional model used is inadequate, under certain operating 
conditions and in regions where abrupt changes in flow area occurs, to describe the three dimensional nature of the flow. 
Another reason is the usage of forced flow correlations for heat transfer coefficient and friction factor in NC flow and 
neglecting the effect of form losses in the various components in the loop. 


Table 3: Maximum absolute errors between experimental and theoretical results for the primary pressure versus NC mass flow rate for 
different power values. 


Max. error (%) 
Power (KW) P 
30 2.1 
60 3.2 
100 0.96 


4. Conclusion 

In this paper, one dimensional continuity, momentum and energy equations are solved to find the mass flow rate of the 
NC flow versus the primary pressure in a PWR plant. It is found that the mass flow rate of the two-phase NC flow is 
decreasing sharply with decreasing the primary pressure. Also, it is found that the NC mass flow rate is slightly dependent 
on core power. The modelling results of this paper are compared with previous experimental data and show a good agreement. 
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It can be concluded that the quasi steady-state one dimensional analyses of this paper can be used to describe the two-phase 
NC phenomena in a PWR during a small break loss of coolant accident. 
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